Imaging charge density fluctuations in graphene using Coulomb blockade spectroscopy 
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Using scanning tunneling microscopy, we have imaged local charge density fluctuations in mono- 
layer graphene. By placing a small gold nanoparticle on the end of the STM tip, a charge sensor is 
created. By raster scanning the tip over the surface and using Coulomb blockade spectroscopy, we 
map the local charge on the graphene. We observe a series of electron and hole doped puddles with 
a characteristic length scale of about 20 nm. Theoretical calculations for the correlation length of 
the puddles based on the number of impurities are in agreement with our measurements. 

PACS numbers: 73.22.Pr, 68.37.Ef 



Graphene is the two dimensional form of carbon con- 
sisting of a single layer of carbon atoms arranged in 
a hexagonal lattice. This simple geometry is a repos- 
itory of exceptional electronic, mechanical, and chemi- 
cal properties. These properties emerge from the lin- 
ear band structure of graphene. The charge carriers 
in graphene behave as massless Dirac fermions and ex- 
hibit very high mobilityii^. This makes graphene very 
promising for device applications^. Efforts in this direc- 
tion have yielded spectacular results from using graphene 
as a nanopore template for DNAf^ to producing 30-inch 
monolayer graphene films on a copper substrate, detach- 
ing them and fabricating transparent electrodes^. How- 
ever, one of the challenges in order to make graphene 
devices is controlling the electrostatic environment. In 
particular, it has been shown that graphene on silicon 
dioxide (Si02) tends to form a series of electron and hole 
puddles near the Dirac point^^. The puddles arise due 
to the random charged impurities on the graphene layer 
and in between the graphene and Si02. This reduces 
the mobility of graphene. The puddles are also partly 
responsible for the minimum conductivity in graphene. 
This limits the on-off ratio for graphene transistors as 
the 'off' state is not reached^ . Hence a complete un- 
derstanding of puddles is vital for progress in graphene 
applications. 

In this work, we use a nanoparticle on the end of a 
scanning tunneling microscope (STM) tip as a sensitive 
charge sensor to spatially map variations in the electron 
charge density. Our technique allows a direct quantita- 
tive measurement of the charge density fluctuations with 
a spatial resolution of several nanometers leading to an 
estimation of the number of impurities. Previous studies 
of electron and hole puddles on graphene lacked the spa- 
tial and/or energy resolution to fully resolve the puddles 
and obtain information about their characteristic distri- 
bution. Charge density variations on a length scale of 
~150 nm have been measured with a single electron tran- 
sistor (SET)^. The size of the SET and distance from the 
graphene limited the spatial resolution of these measure- 
ments. Previous STM measurements of the electron and 
hole puddles lacked the energy resolution of our current 
technique because they relied on a shift in the local den- 



sity of statesii^. In these measurements, we use the sharp 
peaks in dl/dV measurements due to Coulomb blockade 
to significantly improve the energy resolution. 

Graphene was first isolated by exfoliating it from 
3D graphite and attaching it onto a Si02 substrate^. 
There are many alternative ways to make graphene, 
like epitaxy on silicon carbide substrate^ or chemi- 
cal vapor deposition growth on copper— and nickeli^. 
However, the mechanical exfoliation technique still pro- 
vides the highest quality devices. Here, we have stud- 
ied exfoliated graphene on a Si02 substrate using a 
STM operating in ultra-high vacuum at 4.6 K. Previ- 
ous STM studies on exfoliated graphene have been in- 
strumental in correlating its structure with the electronic 
propertiesii^ii^. In this work, we report density fluctua- 
tions in graphene using Coulomb blockade spectroscopy 
measurements. Coulomb blockade is the suppression 
of electron tunneling when the charging energy e^/2C, 
where C is the capacitance of the nanoparticle, is unavail- 
able to the electrons at low temperature and voltages, i.e. 
ksT, eV « e^/2Gi^. With a nanoparticle at the end of 
the STM tip, the system exhibits the two junctions which 
are necessary to produce Coulomb blockade. The tun- 
neling gap between the nanoparticle and the surface acts 
as one junction, while the barrier between the tip and 
nanoparticle is the other. The high spatial resolution of 
the STM and the low temperature allows a spatially re- 
solved measurement of Coulomb blockade. As Coulomb 
blockade is sensitive to not only changes in capacitance 
but also the electrostatic environment of the nanoparti- 
cle, we are able to detect charge fluctuations on graphene 
using Coulomb blockade spectroscopy. 

Graphene flakes were exfoliated from graphite on a Si 
substrate with 300 nm of thermally grown oxide. Ti/Au 
electrodes were then deposited on the graphene through a 
shadow mask^^. This technique does not require the use 
of PMMA as a resist. Therefore it yields clean graphene 
devices that need no additional cleaning procedures prior 
to imaging with STM. The graphene sample is trans- 
ferred to the ultra-high vacuum (p < 10~^^ mbar) STM 
and cooled to 4.6 K. Electrochemically etched tungsten 
tips were used for imaging and spectroscopy. A constant 
density of states for the tip was confirmed on a Au sur- 
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FIG. 1. (color online) (a) Schematic of the graphene sam- 
ple with gold contacts on a Si02 substrate. The nanoparticle 
at the tip apex is shown in the close up sketch, (b) Optical 
microscope image of the monolayer graphene flake and gold 
electrodes, (c) Large scale STM image of the graphene sur- 
face. (0.5 V, 100 pA, 100 nm x 100 nm). The scale bar is 
20 nm. (c) STM image showing the atomic resolution of the 
hexagonal lattice of monolayer graphene (0.5 V, 100 pA, 10 
nm X 10 nm). The scale bar is 2 nm. 

face. After confirmation of the density of states of the 
tip, a nanoparticle was attached to the end of it and spec- 
troscopy measurements on graphene were performed. 

Figure [TJa) shows a schematic diagram of the mea- 
surement setup. The bias voltage is applied to the STM 
tip and the sample is grounded. The particle on the tip 
apex is a Au nanoparticle attached after the tip condi- 
tioning on the Au surface. Because of the cleanliness of 
the samples, we are able to obtain large area images of 
the graphene surface which are free of defects or residue 
from sample processing. Figure [Hb) shows the optical 
microscope image of the monolayer graphene flake with 
gold contacts on the Si02 substrate. Figure [TJc) shows 
a typical 100 x 100 nm area of the graphene surface. 
The sample is smooth with a slowly varying height due 
to the underlying Si02 substrate. There are no sharp 
features that arise from the presence of contaminants on 
the surface or defects. Atomic resolution topography of 
the sample is shown in Fig. md) showing the hexago- 
nal lattice of carbon atoms as expected for monolayer 
graphene along with height variation due to the uneven 
Si02 substratei^^^^. 

We have measured the local density of states (LDOS) 
as a function of energy. For this measurement, the tip is 
held at a fixed location and the feedback loop is turned 
off. The voltage between the tip and substrate is ramped 
and the differential conductance as a function of voltage 
is acquired. An ac modulation voltage of 5 mV rms at 
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FIG. 2. (color online) (a) dl/dV spectroscopy showing peaks 
due to the addition of individual electrons to the nanoparti- 
cle. (b) dl/dV as a function of tip voltage and gate voltage, 
showing a series of diamonds. This demonstrates that the 
nanoparticle on the end of the tip acts as a quantum dot. (c) 
dl/dV curves taken at two different points on the graphene 
surface. Each shows a series of peaks but they are horizontally 
offset due to the different potential on the graphene surface, 
(d) Distribution of peak spacings for a 100 nm x 100 nm 
area of the graphene film. The peaks are clustered around 
the same value showing that there is a single quantum dot 
located on the end of the tip. 



570 Hz is applied to the tip and the resulting current is 
detected with a lockin amplifier. Figure [2] (a) shows a 
typical dl/dV curve acquired with a nanoparticle at the 
end of the tip. There are a series of sharp peaks that are 
equally spaced in energy. This is in contrast to a mea- 
surement on graphene without the nanoparticle which 
shows a smoothly varying density of state^ii^. The 
peaks arise due to the addition of individual electrons 
to the nanoparticle. The setpoint current and voltage 
sets the overall resistance between the tip and graphene. 
Using a small setpoint current ensures that the resistance 
between the nanoparticle and graphene, Ksub^ is always 
much greater than the resistance between the nanopar- 
ticle and tip, Kup- Because of this strong asymmetry in 
the tunnel barriers, the barrier between the nanoparticle 
and graphene is always the rate limiting barrier. There- 
fore, we obtain a Coulomb staircase with the spacing of 
the peaks determined by the capacitance between the 
nanoparticle and graphene^^. 

We have measured the dl/dV curves for different back 
gate voltages on the monolayer graphene sample. Once 
again the curves are characterized by sharp peaks which 
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vary as a function of gate voltage. Figure [2] (b) plots the 
dl/dV curves as a function of gate voltage and tip volt- 
age. It shows the characteristics of a Coulomb blockade 
diamond plot providing further evidence of the presence 
of a small nanoparticle on the tip. In the region around 
zero sample voltage, the bright red diamond shaped re- 
gions show that there is no flow of current, indicating 
Coulomb blockade. As the voltage increases a series of 
peaks appear. Between each peak the number of elec- 
trons on the nanoparticle changes by one. The gate acts 
to induce charge on the nanoparticle and therefore the 
peaks change in energy. Only one set of peaks is strong, 
this is due to the asymmetry of the tunnel barriers. This 
strong peak corresponds to lining up the energy of a state 
on the nanoparticle with the energy of the graphene. 
Therefore, as the voltage between the tip and graphene 
is increased, we see a series of peaks separated by e/Csuh 
where Csuh is the capacitance between the nanoparticle 
and graphene. The strong peaks slope towards decreas- 
ing tip voltage as the gate voltage is increased. This is the 
opposite direction expected for a normal gate electrode 
which would tend to induce negative charges as the gate 
voltage increases. However, the effect of the gate on the 
nanoparticle is screened by the graphene. A positive volt- 
age on the gate induces electrons on the graphene which 
then tend to shift the energy levels on the nanoparticle. 
Effectively, the graphene layer acts to reverse the sign of 
the gate voltage. 

As the tip is moved to different positions of the 
graphene, the energies where peaks occur shift. This 
is seen in Fig. [2] (c) which shows the spectroscopy re- 
sults from two locations on the graphene surface sepa- 
rated by a distance of 10 nm. The entire dl/dV curve 
has been shifted along the tip voltage axis. All of the 
peaks still have the same spacing showing that the capac- 
itance between the nanoparticle and graphene remains 
unchanged. However, the energy at which the peaks oc- 
cur has changed. This is due to the presence of an offset 
charge on the nanoparticle. Below we will show that this 
offset arises due to variations in the potential on the sur- 
face of the graphene. Figure [2] (d) shows the distribution 
of peak spacings over a 100 nm x 100 nm area. As the 
peak spacings cluster around a single value, this shows 
that the we are using the same nanoparticle for all of our 
measurements. On the other hand, the distribution of 
the energies of the peaks over the same area shows a uni- 
form distribution. When a new nanoparticle is attached 
to the tip, the peak spacing changes but otherwise the 
results are unchanged. 

Figure [3] shows the local density of states over a 100 
nm by 100 nm area of the sample. Six images from a 
movie of the density of states as a function of tip voltage 
are shown^^. The bright regions represent the locations 
where electrons are added to the nanoparticle. They are 
peaks in the dl/dV curve. Following one of the bright 
lines corresponds to moving along a path of constant in- 
duced charge on the nanoparticle. As the voltage be- 
tween the tip and graphene is changed, the location where 
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FIG. 3. (color online) Local density of states maps taken at 
a series of tip voltages as indicated for each map. Each map 
is 100 nm x 100 nm. A ring is circled to show how it changes 
with energy. They start out circular (a-b) and slowly become 
asymmetric (c-e). Then a second ring forms (f) at the same 
location once the charge on the nanoparticle changes by one. 



electrons are added to the nanoparticle also changes. Fig- 
ure [3] (a) to (f) shows the evolution of the rings as the tip 
voltage is changed. Notice that the circled ring grows as 
the voltage goes from 0.47 V to 0.37. When the sample 
voltage has reached 0.30 V, Fig. [3](f), an additional elec- 
tron has been added to the nanoparticle and the image is 
equivalent to Fig. [3] (a) at 0.47 V. The spatial variation in 
the density of states map implies that the electrostatic 
environment of the nanoparticle is changing as a func- 
tion of position. Therefore, these maps give information 
about the local potential or charge density fluctuations 
in graphene as discussed below. 

To obtain a map of the local potential or equivalently 
the charge density, we can analyze the energy at which 
the peaks occur as a function of position. Since the peaks 
are arising due to Coulomb blockade, we expect that the 
spacing in energy between peaks will always remain con- 
stant at e/Csub' This is the reason that Fig. [3] (a) and 
(f) look identical; they are separated in energy by e/Csub- 
On the other hand, the energy at which the peaks occur 
gives information about the amount of charge which must 
be induced on the nanoparticle in order for it to conduct. 
If the offset charge on the nanoparticle is e/2, then we 
expect a series of peaks at ne / Csub where n is the number 
of electrons induced on the nanoparticle. As the offset 
charge changes, the energies of the peaks will shift such 
that they occur at ( AQ + ne) / Csub where AQ is the off- 
set charge relative to e/2. In both cases, the peaks are 
spaced by e/Csub- We have used the data in Fig. [3] to 
measure the offset charge at each point on the graphene 
surface. This offset charge can be converted to a potential 
difference, AV, by dividing by Csub- Furthermore, the 
potential difference can be converted to a local electron 
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FIG. 4. (color online) Charge density fluctuations in the same 
100 nm by 100 nm region of the flake as in Fig. [S] Red areas 
are positively charged while blue areas are negatively charged. 

density using the band structure of graphene. Using the 
linear dispersion relation of graphene, we find that fluc- 
tuations in charge density are related to fluctuations in 
potential by An = {/SVY /irihvFY where vp = 1.1 x 10^ 
m/s is the Fermi velocity. The results are shown in Fig. 
m There are regions of positively and negatively charged 
puddles on the surface of the graphene. This map shows 
a typical size scale of electron and hole puddles of 20 
nm. The regions of positive charge correspond to areas 
of the graphene flake where the rings in Fig. [3] are shrink- 
ing as the tip voltage is decreased. Likewise, rings that 
grow with decreasing tip voltage, such as the ring in the 
pink ellipse in Fig. [3] are areas of negative charge on the 
graphene. 

We find that the FWHM of our distribution of charge 
densities is An = 3 x 10^^ cm~^. This is about a factor 
of 6 higher than found in previous SET measurements 
of electron and hole puddles on graphene^. However, 
our significantly improved spatial resolution leads to the 
higher variation because we are now sensitive to charge 
fluctuations on the nanometer scale as opposed to the 
previous measurements which were limited to about 150 
nm resolution. With the lower spatial resolution, the 
charge fluctuations tend to be averaged away. Our value 
for the charge fluctuation is similar to previous STM mea- 
surements of electron and hole puddles^, however this 
technique has the ability to resolve much smaller changes 
in charge density because of the high energy sensitivity 
of the Coulomb blockade spectroscopy. 

By performing an auto-correlation of the image in Fig. 
m we can find the correlation length of the puddles. The 
results are shown in Fig. [5] (a). It shows a FWHM of 
the distribution of 20 nm. Theoretical calculations have 
shown that scattering from random charged impurities 
causes electron and hole puddles and their correlation 



length varies with the concentration of these impuritie^. 
We can estimate the number of impurities in our sample 
from the electrical conductivity at high gate voltages^!. 
The conductivity is given by a = 2{)eeVg / {hriit) where 
Ui is the impurity density, e is the dielectric constant of 
Si02, t is the oxide thickness and Vg is the gate voltage. 
Figure [5] (b) shows the resistivity as a function of gate 
voltage. From this plot, we can estimate the density of 
impurities as n^ ^ 1 x lO^^cm"^. This impurity concen- 
tration gives a theoretical correlation length for the pud- 
dles of 12 nm. The global nature of the transport mea- 
surement compared to our local measure of the charge 
density fluctuations may explain the small discrepancy 
between the theoretical puddle size and our images. 

Thus we have shown that a nanoparticle at the STM 
tip apex can act as a sensitive charge sensor to detect the 
charge density fluctuations in graphene. Using spatially 
resolved spectroscopy measurements we have mapped 
these fluctuations and found a characteristic puddle size 
of 20 nm. A comparison with theoretical calculations 
shows that our graphene samples have an impurity den- 
sity on the order of 10^^ cm~^ which give rise to these 
puddles. In order for graphene devices to reach their full 
potential the number of impurities must be significantly 
reduced. This technique of using a gold nanoparticle at 
the tip apex for spatially resolved measurements will en- 
able us to detect the reduction in the number of impu- 
rities and consequently test the quality of the flakes for 
device fabrication. 
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FIG. 5. (a) Correlation length for the puddles shown in Fig. 
[41(b) Resistivity as a function of gate voltage for the graphene 
flake. 
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